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ABSTRACT 

! We report the first detection of a planetary-mass companion to a star using the High Resolution 

Spectrograph (HRS) of the Hobby-Eberly Telescope (HET). The HET-HRS now gives routine radial 
velocity precision of 2-3 ms~^ for high SNR observations of quiescent stars. The planetary-mass 
companion to the metal-rich KOV star HD 37605 has an orbital period of 54.23 days, an orbital 
eccentricity of 0.737, and a minimum mass of 2.84 Jupiter masses. The queue-scheduled operation 
of the Hobby-Eberly Telescope enabled us to discovery of this relatively short-period planet with a 
total observation time span of just two orbital periods. The ability of queue-scheduled large-aperture 
^— ^ ' telescopes to respond quickly to interesting and important results demonstrates the power of this 

' new approach in searching for extra-solar planets as well as in other ares of research requiring rapid 

, response time critical observations. 

Subject headings: stars: individual(HD 37605) — techniques: radial velocities — planetary systems 

Q ■ 1. INTRODUCTION Hobby-Eberly Telescope, (HET) permits the observer to 

ti '. Traditional ground based radial velocity searches for '"eact quickly to the data and to obtain timely foUow- 

^ ■ extra-solar planetary systems have used medium- and observations of critical targets at the orbital phases 

^ ■ large-size telescopes assigned to the program for particu- ^^^re they are needed. We report use of the HET to 

lar nights of observing. While this type of approach has 'i^tf * ^ P "f^^ oZ^?™Sf'''°u.','^ * t , .^r.l''''^ 

been extremely successful, resuhing in about 120 detec- ^^^^ star HD 37605. The abihty of the HET observing 

, tions of planetary-mass companions to solar-type stars, ^^^^^^ accommodate our observing needs enabled us to 

H ■ the vagaries of the telescope scheduling process intro- discover the planet with only a total of 100 days (slightly 

■ - ■ duces certain restrictions and constraints on the opti- ^^ss than two orbital periods) elapsing between the first 

mal acquisition of time-critical observations. In radial ^,1^^ ^^^^1 observations. We discuss both the detection of 

velocity surveys, it is often the case that once a candi- ^he planetary companion itself and the special attributes 

date planetary system has been identified and a prelim- ,the HET that enabled the fina orbital confirmation 

inary orbit determined, additional observations must be ^e obtained m a very timely and expeditious manner, 
obtained at particular orbital phases in order to place 



> 



■ 

X 



, . 1 , 4- ■ 4- 1 1 1 r^u 2. THE HOBBY-EBERLY TELESCOPE AND ITS HIGH 

tighter constraints on particular orbital elements. Ob- RESOLUTION SPECTROGRAPH 
taining the necessary data at the optimal orbital phase 

may be extremely difficult if observing runs are scheduled The Hobby-Eberly Telescope consists of an array of 91 

for a few nights once per month, or even less frequently. hexagonal shape, spherical figure mirrors, each of 1 meter 

Thus a telescope that is queue scheduled, such as the diameter,^ in a fixed truss tilted at 35° zenith distance, 

resulting in a maximum effective aperture of 9.2 m diam- 

*BASED ON OBSERVATIONS OBTAINED WITH THE eter. Declination ranges in the sky between -11° and 

HOBBY-EBERLY TELESCOPE, WHICH IS A JOINT +72° ^ay be selected by rotating the telescope truss in 
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UNIVERSITY, LUDWIG MAXIMILIANS UNIVERSITAT bundle along the telescope focal plane while the tele- 

MtJNCHEN, AND GEORG AUGUST UNIVERSITAT scope structure remains fixed, and feeding the light into 

GOTTINGEN. the spectrograph. The HET design was optimized for 
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employing different instrument configurations. The HET 
data distribution pohcy is to allow the PI to gain ac- 
cess to the pervious nights data every morning allowing 
rapid feed back to the HET astronomers. This is an ideal 
mode of operation for our high-precision radial velocity 
program. 

The High Resolution Spectrograph (HRS) ( TuU 199|) 
for the HET was built by a team headed by Robert Tull 
and Phillip MacQueen. This spectrograph was designed 
to be able to make stellar radial velocity variation mea- 
surements with a precision of 3ms~^ or better on stars 
as faint as V — 10. The instrument is housed in an in- 
sulated enclosure in the basement of the HET building, 
and thus is not moving with the telescope. The image 
scrambling provided by the 34.3 m long optical fiber feed 
helps ensure that the pupil seen by the spectrograph is 
always illuminated in a consistent manner. Only a set 
number of instrument configurations are available, and 
these are all defined by kinematic mounts for the grating 
and cross-disperser positions. The I2 gas cell is contained 
in a coUimated section in front of the spectrograph slit, 
so refocusing of any optics is unnecessary for I2 cell use. 
The detector is a mosaic of two 2048 x 4100 pixel E2V 
CCDs. For our high-precision radial velocity observa- 
tions, we use the HRS in a mode giving resolving power 
of 60,000. The echelle grating was used on the center 
of the blaze, and the cross disperser was positioned so 
that the break between the "red" and the "blue" CCD 
chips was at around 5940A. For most seeing conditions, 
a 2 arcsec diameter optical fiber was used. 

So far, we have observed 173 F-M dwarfs four or more 
times. The F through K stars are selected to be stars 
with low levels of stellar activity based on X-ray emission, 
photometric variability, and measured Ca II H&K emis- 
sion. We intentionally do not bias the selection on stel- 
lar metallicity, but we do attempt to exclude stars with 
usini greater than about 15kms~^ and known short- 
period binary stars. The separate M-dwarf survey an d 
its selection criteria are discussed bv lEndl "eTall ll200l . 
Of all of these stars, 20 show large-amplitude variations 
indicative of previously unknown binary star systems, 
and 11 additional stars show rms variations greater than 
20ms~^, but probably not large enough to be due to bi- 
nary stellar companions. These stars represent good can- 
didates for short-period planetary companions. Figure ^ 
shows a histogram of the velocity rms about the mean 
of our time-series measurements for the remaining 142 
stars having an rms of 20ms~^ or less, after any statisti- 
cally significant linear velocity trends have been removed. 
The time span of the data for each star ranges from one 
week to greater than three years. The median rms of 
the 123 stars in Figure ^ with rms of 13 ms""'^ or less is 
4.0ms^^. If we examine just the 112 stars with rms of 
8ms~^ or less (the bulk of the sample), then both the 
mean and the median rms of the remaining "stable" stars 
are just 3.6ms^^. The peak of the histogram is centered 
in the 2-3 m s^^ bin. We note that a significant contribut- 
ing factor to the observed velocity rms is intrinsic to the 
stars themselves. We have not attempted to correct these 
rms values in any way for velocity variations induced 
by stellar photospheric and chromospheric activity ("jit- 
ter"). In general, the stars for which we have obtained 
the lowest rms are old, inactive stars with spectra of 
high signal-to-noise ratio. Thus, Figure ^ clearly demon- 
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Fig. 1. — Histogram of the velocity rms of 142 stars observed 
in this program with the Hobby-Eberly Telescope High Resolution 
Spectrograph. 



strates that we have achieved routine radial velocity pre- 
cision of 3ms~^ or better with the HET-HRS in routine 
"production mode" on a large sample of a wide variety 
of stars over a substantial time period. Indeed, a sig- 
nificant fraction of our stars are showing radial velocity 
precision of 2 ms~^ or better, which is comparable to the 
precision achieved f or high SNR observ ations with the 
ESO VLT /UVES bv lKilTster et al.l ll2003^. Our2-3ms-i 
precision is significantly bet ter than that obta ined with 
COR ALIE ijNaef et al.i20 011 or ELODIE ( Bar anne et alJ 
Il996j) . and is directly competitive with the 3 ms~^ pre- 
cision which i s achieved on the Keck HIRES l|Vogt et alJ 
I2000|) . AAT i Jones et al.ll20 02^. or the 2ms-i demon- 
strated by HARPS llPepeet a l. 2004). 

3. THE PLANETARY COMPANION TO HD 37605 

3.1. Characteristics of the host star 

HD 37605 (HIP 26664, BD-l-05 985, SAO 1130 15, G 
99-22 ) is a V=8.69 KOV star. The Hipparcos ijESAI 
I1997|) parallax of 23.32 mas corresponds to a distance of 
42.9parsec, and gives an absolute magnitude of My = 
5.51. The stellar parameters: effective temperature, Teff, 
surface gravity, logf/, and microtur bulence, £ list ed in 
Tabled were derived using MOOG IISnedenlll97.'f with 
atmosp heres based on the 1995 version of the ATLAS9 
code (Castelli et al. 1997) following the procedure out- 
lined in Pa ulson et al.. (2003) and briefly described be- 
low. Within IRAF we fit the 20 Fel an d 12 Fell lines 
listed in Table 1 of iPaulson et al.l l|2003f) with Gaussian 
profiles to obtain equivalent widths. The T^g was de- 
rived by requiring that individual line abundances be in- 
dependent of excitation potential. Likewise, ^ was de- 
rived by requiring that individual line abundances be in- 
dependent of line strength. The logg was determined 
by forcing ionization equilibrium between neutral and 
singl y ionized Fe. We adopted solar loge(Fe) of 7.52 
from lAnders fc Greve ssc (198^, which yields an [Fc/H] 
for HD 37605 of 0.39±0.06. 
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TABLE 1 
Stellar properties for, HD 37605 



200 



Parameter 




derived value 


Tcff 


5475± 


50K 


logs 


4.55± 


0.1 




0.8± 


0.2 km s-l 


[Fe/H] 


0.39± 


0.07 



3.2. HET Radial Velocity Observations 

Our typical HET observing technique is to obtain 
about 4-5 initial radial velocity measurements over the 
course of 1-2 weeks, in order to search for short-period 
"hot- Jupiter" RV variability. If a star appears stable on 
short timescales, it is then scheduled for less frequent 
observations in order to search for longer period orbits. 
For HD 37605 {V = 8.7), exposure times were 900 sec- 
onds, giving a typical signal-to-noise ratio of about 250 
per resolution element. The initial set of observations 
of HD 37605 were constant to within the observational 
error, but the next observation taken about one month 
later showed a decrease of about 200 ms"^. We then 
took advantage of the queue scheduled operations of the 
HET to put the star back into the queue at high prior- 
ity for frequent observations. Following our request for 
frequent coverage, the HET Resident Astronomers then 
obtained spectra every 3-4 days during the decrease in 
radial velocity. Data from each night's observations were 
reduced and velocities computed the following morning. 
We were easily able to fit a new orbital solution as the 
data from each night became available. It quickly be- 
came obvious that the radial velocity minimum and pe- 
riastron passage would occur just as the star was being 
lost from the HET observability window; the last HET 
data points would be obtained during evening twilight. 
It was critical to attempt to obtain nightly HET veloc- 
ity measurements during this orbital phase, as the depth 
of the velocity minimum would constrain the orbital ec- 
centricity and the K velocity (and hence the companion 
minimum mass). Since these were going to be difficult 
HET observations, and it was quite possible that the 
HET would no longer be able to observe the star be- 
fore the periastron passage was complete, we arranged 
for the McDonald Observatory 2.1m telescope observer 
to obtain radial velocity measurements with the Sandi- 
ford Cassegrain Echelle spectrograph ijMcCarthv et al.| 
11993^ and its I2 cell during this interval. This instrument 
records the 5000 to 6OOOA region at resolving power of 
60,000. These 2.1m data would be of significantly lower 
velocity precision due to the much lower signal-to-noise 
ratio and the inherent mechanical and thermal instability 
of a large heavy Cassegrain spectrograph, but they would 
certainly be adequate to cover this crucial interval when 
the HET might not be able to observe the star. The HET 
data for HD 37605 are shown as filled dots in Figure El 
and the 2.1m Sandiford echelle data are shown as open 
triangles. The observed velocities from both telescopes 
are given in Table [3 The velocities in this table have 
been corrected for the different (and arbitrary) velocity 
zero points of the two instruments used. The simulta- 
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Fig. 2.— Radial velocity data for HD 37605 from the HET HRS 
(filled circles) and from the McDonald Observatory 2.1m Sandiford 
Cassegrain echelle (open triangles). The size of most of the HET 
error bars is smaller than the size of the data points. The best-fit 
combined Keplerian orbital solution is plotted as solid line. The 
rms residual of the HET/HRS data from the orbital solution is 
4.7ms~^. 



TABLE 2 

Differential radial velocities for HD 37605 



JD-2400000 


velocity 


uncertainty 


instrument 




(m s"-') 


(m s"-') 




53002.6715 


94.7 


8.4 


HET HRS 


53003.6853 


97.1 


7.8 


HET HRS 


53006.6620 


92.8 


7.0 


HET HRS 


53008.6641 


98.6 


6.2 


HET HRS 


53010.8048 


90.8 


6.2 


HET HRS 


53013.7940 


79.3 


6.6 


HET HRS 


53042.7280 


-126.4 


7.0 


HET HRS 


53061.6676 


101.7 


4.2 


HET HRS 


53065.6468 


94.8 


4.8 


HET HRS 


53071.6438 


67.5 


4.6 


HET HRS 


53073.6382 


64.9 


4.2 


HET HRS 


53082.6237 


24.6 


4.1 


HET HRS 


53083.5954 


18.2 


5.6 


HET HRS 


53089.5958 


-19.0 


4.0 


HET HRS 


53092.5980 


-47.9 


4.1 


HET HRS 


53094.5866 


-71.2 


3.3 


HET HRS 


53095.5864 


-89.6 


3.5 


HET HRS 


53096.5874 


-107.7 


6.7 


HET HRS 


53098.5763 


-190.0 


4.2 


HET HRS 


53102.5727 


-409.9 


5.5 


HET HRS 


53072.1363 


69.5 


25.99 


2.1m 


53097.6345 


-177.0 


30.24 


2.1m 


53098.6555 


-181.7 


27.91 


2.1m 


53101.6647 


-412.1 


78.12 


2.1m 


53102.6034 


-297.5 


23.58 


2.1m 


53103.5928 


-203.0 


20.37 


2.1m 


53104.5972 


-79.7 


29.20 


2.1m 



neous orbital solution to both data sets is also shown 
in Figure El The rms residual of the HET/HRS data 
from this orbital solution is 4.7 ms~^. The rms residual 
of the 2.1m Sandiford echelle data from the solution is 
44.2ms~^ 

3.3. Orbital Solution 

We performed a simultaneous fit of the HET-HRS and 
2.1m Sandiford echelle velocity data using the Gauss fit 
ijMcArthur et alJll994^ software, which uses a robust es- 
timation method to find the combined orbital solution, 
where the arbitrary velocity zero points of both data sets 
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TABLE 3 

Orbital parameters for, the companion to HD 37605. 



Parameter derived value 



p 


days] 


54.23± 


0.23 


T 


JD] 


2452994.27± 


0.45 


K 


ms-i] 


262. 9± 


5.5 


e 




0.737± 


0.010 


UJ 


degrees] 


211. 6± 


1.7 


f{m) 


solar masses] 


3.13 X 10-»± 


5.9 X 10""' 


msin j 


Mjup] 


2.84± 


0.26 


a 


AU] 


0.26± 


0.01 



are included as fit parameters. The combined orbital so- 
lution is given in Table 13 and is shown as the solid line 
in Figure 121 The planet is in a highly eccentric orbit. Its 
eccentricity of 0.737 is exceeded only by HD 80606b (e 
= 0.93) and HD 222582b (e=0.76). The periastron dis- 
tance of 0.070 AU is large enough for the planet easily to 
have avoided tidal circularization over the lifetime of the 
system. 

The shape of the observed radial velocity curve is very 
difficult to explain by any mechanism other than orbital 
motion. Nevertheless, we have examined the available 
data to ensure that our interpretation of the observed 
velocity variations as resulting from orbital motion is in- 
deed correct. The Hipparcos photometry of HD 37605 
shows an rms variation of only 0.0020 magnitudes, which 
is comparable to the Hipparcos photometric precision. 
Thus, the observed RV variations can not be linked to 
any photometric variability of the primary star. We have 
also searched for evidence of chromospheric variability of 
HD 37605. Normally, we would prefer to measure the 
Ca H S index from the Ca H H and K lines, but the 
HET primary mirror coatings and the HET optical fibers 
result in a short-wavelength cutoff of the HET HRS at 
around 4200A. We do have one spectrum of HD 37605 ob- 
tained with the McDonald Observatory 2.7m telescope 
2dcoude spectrometer, which includes the Ca H H & K 
lines . From this spe ctrum we measure a McDonald S in- 
dex ijPaulson et alii2002) of 0.158 ± 0.022. This 
comparable to the value found in the old, inactive star 
T Ceti (SmcD = 0.166 ± 0.008), thus indicating a low 
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level of chromospheric activity for HD 37605. We have 
also examined the profiles of the H a absorption line for 
any possible variability that might be related to stellar 
activity. A careful alignment and normalization of the 
H a line profile from all of our regular program spectra 
of HD 37605 shows no variations beyond what is expected 
from photon noise. Thus we conclude that the observed 
radial velocity variations are indeed due to a planetary- 
mass companion moving in a highly eccentric orbit. 

4. CONCLUSIONS 

We report the first discovery of an extra-solar planet 
using the Hobby-Eberly Telescope and its High Res- 
olution Spectrograph. The planetary-mass companion 
to the metal-rich star HD 37605 is in a highly eccen- 
tric orbit. We have demonstrated that the HET High 
Resolution Spectrograph can produce a long-term ve- 
locity precision of 3ms~^ or better over a multi-year 
timespan. The queue-scheduled operation of the Hobby- 
Eberly Telescope made the detection and characteriza- 
tion of the orbit far simpler than it would have been 
with a conventionally scheduled telescope. Once the pos- 
sible orbital motion was detected, we were able to obtain 
immediate follow-up observations at the particular or- 
bital phases that were crucial for the orbit determination. 
This highly flexible mode of telescope operation, coupled 
with the large telescope aperture and the extremely sta- 
ble High Resolution Spectrograph, makes the HET- HRS 
an ideal and extremely efficient and flexible facility for 
large surveys for extrasolar planetary systems. 



The Hobby-Eberly Telescope (HET) is a joint project 
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